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Pb(Feo. 5 Nbo. 5)03 (PFN), one of the few relaxor multiferroic systems, has a G-type antiferromag¬ 
netic transition at Tn = 143 K and a ferroelectric transition at Tc = 385 K. By using high-resolution 
neutron-diffraction experiments and a total scattering technique, we paint a comprehensive picture 
of the long- and short-range structures of PFN: (i) a clear sign of short-range structural correlation 
above Tc, (ii) no sign of the negative thermal expansion behavior reported in a previous study, and 
(hi) clearest evidence thus far of magnetoelectric coupling below Tn. We conclude that at the heart 
of the unusual relaxor multiferroic behavior lies the disorder between Fe^^ and Nb^^ atoms. We 
argue that this disorder gives rise to short-range structural correlations arising from O disorder in 
addition to Pb displacement. 

PACS numbers: 75.85.+t,61.05.F-,75.40.-s 


INTRODUCTION 

It is rare in nature to find a single system which hosts 
two or more ordered ground states out of otherwise un¬ 
connected degrees of freedom. If multiple ordered ground 
states exist, a natural question to ask is whether they are 
connected to one another. Although the same electrons 
may be responsible for multiple degrees of freedom, allow¬ 
ing coupling of the various forms of order, this is uncom¬ 
mon, making materials in which they are coupled particu¬ 
larly interesting. Magnetoelectric multiferroic materials, 
in which ferroelectricity and magnetism coexist and are 
coupled, are an example of just such a material. The 
ability to control one form of order via another offers a 
huge potential for technological applications and, at the 
same time, poses new challenges for our understanding 
of how distinct degrees of freedom as disparate as bulk 
polarization and magnetization can be coupled to one 
another IDE]. Here the sought-after coupling, so-called 
magnetoelectric effects, can lead to better manipulation 
of unusual multiferroic behavior and to exotic excitations 
of mixed character. Despite the huge interest, however, 
the origin of magnetoelectric coupling has often proven 
challenging to address experimentally for a given mate¬ 
rial. 

Lead iron niobate Pb(Feo. 5 Nbo. 5)03 (PFN) is a 
multiferroic material with a ferroelectric transition at 
Tc = 385 K [3] and an antiferromagnetic transition at 
Tn = 143 K [3]-[6] which is known to be G-type [7]. It is 
noteworthy for its high dielectric constant, which changes 
at both the ferroelectric man] and the magnetic tran¬ 
sitions laigiiQ] and is strongly frequency dependent. Its 
reported high dielectric constant makes it a suitable can¬ 


didate material for multilayer ceramic capacitors among 
other electronic devices. The reported strong frequency 
dependence, making it a rare example of a relaxor multi¬ 
ferroic, most likely arises from disorder at the magnetic B 
site of the perovskite structure. Since many ferroelectric 
systems have disordered magnetic ions on the B site, a 
full understanding of PFN might well have wider implica¬ 
tions for finding or optimizing other relaxor multiferroic 
materials for potential applications. 

Despite the interest in the physical properties of PFN, 
however, several features of the underlying crystal struc¬ 
ture still remain unresolved. For example, two compet¬ 
ing space groups, R3m [7] and Cm [niiia, have been 
proposed for the low-temperature structure, and there is 
an unconfirmed report of negative thermal expansion be¬ 
low the antiferromagnetic transition m- In attempting 
to characterize the relaxor behavior, several groups have 
investigated the local structure [niiTij. which is closely 
associated with relaxor behavior in ferroelectrics, and dis¬ 
cussed the possibility of structural disorder both exper¬ 
imentally [HI US] and theoretically m- These studies 
notwithstanding, the key questions, we believe, are yet 
to be fully answered. For instance, it will be interest¬ 
ing to know how the relaxor behavior is connected to 
the short-range local structure. More importantly, one 
first has to know the exact details of the local distortion 
and, if possible, the structural basis of the sought-after 
magnetoelectric coupling. 

In order to address these questions, we have un¬ 
dertaken full high-resolution neutron powder-diffraction 
studies as well as total scattering experiments using two 
state-of-the-art instruments: S-HRPD and NOVA, both 
at J-PARC, in Tokai, Japan. For this kind of study it is 
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essential to be able to examine both long- and short-range 
structures simultaneously as performed for other ferro¬ 
electric materials m- By combining data obtained from 
both instruments, we conclude that the low-temperature 
space group is Cm with clear signs of the short-range 
structure surviving even above the ferroelectric transi¬ 
tion temperature, offering an explanation for the relaxor 
behavior. Neither experiment showed signs of negative 
thermal expansion. We then uncovered structural signa¬ 
tures of the magnetoelectric coupling by analyzing the 
temperature dependence of the electric polarization cal¬ 
culated from the structure parameters. Our detailed 
analysis of the local structure leads us to the conclusion 
that both Pb displacement and O disorder are exhibited 
in this material. 


EXPERIMENTAL DETAILS 

About 7 g of polycrystalline PFN samples were pre¬ 
pared by a standard solid-state method. Stoichiometric 
PbO, Fe 203 , and Nb 205 were mixed, then calcined in 
air at 850 °C for 2 days. After calcination the products 
were ground, pressed into pellets, and then sintered at 
1050 °C for 1 day. The samples were verified to be single 
phase with a Rigaku Miniflex II x-ray diffractometer as 
well as the later high-resolution neutron diffraction stud¬ 
ies which form the main body of this paper. 

High-resolution neutron time-of-fiight powder- 
diffraction experiments were carried out using the 
S-HRPD beamline at J-PARC, Tokai, Japan, on a 
powder sample in a cylindrical vanadium can, at ten 
temperatures from 10 to 300 K. Rietveld refinement of 
the diffraction data was performed using FULLPROF [18] . 
The subsequent total scattering experiments were 
carried out using a vanadium-nickel sample can at the 
NOVA beamline at J-PARC with a maximum Q value 
of Qmax = 30 at five temperatures between 62 and 
453 K. Sample environment constraints prevented mea¬ 
surements above the ferroelectric transition temperature 
at S-HRPD, so NOVA was also used to collect diffraction 
patterns at selected higher temperatures. For the local 
structure study, a radial distribution function was 
calculated from the total scattering data and analyzed 
using PDFGUI [19]. The experimentally obtained S{Q) 
was then Fourier transformed into real space to obtain 
the atomic pair distribution function G(r) as follows: 

O rQmax 

G{r) = 47rr [p{r) - po] = - / [^'(Q) - 1] sinQrdQ, 

^ Jo 

where p(r) and po are the atomic number density and 
the average number density, respectively, and Qmax is 
the maximum scattering wave vector. G{r) was modeled 
using 4x4x4 supercells for several different models of 
the short-range structure. 


RESULTS AND DISCUSSION 



FIG. 1. (a) Neutron powder-diffraction patterns with refined 
results in the Cm and R3m (inset) space groups, respectively, 
at 300 K. Circles represent data, a line the refinement, and the 
line at the bottom the residual, whereas green vertical bars 
mark the Bragg positions, (b) Diffraction patterns at several 
temperatures. The upper two datasets (starred) are NOVA 
results, and the others are S-HRPD results. Several peaks 
disappear above Tc = 385 K due to the structural transition 
(arrows). 

The low-temperature space group of PFN is thought 
to be the monoclinic space group Cm, but as discussed, 
the R3m space group has also been proposed as an al¬ 
ternative structure; to the best of our knowledge it is 
still not settled which one of the two space groups is cor¬ 
rect. Figure [^a) shows the neutron powder-diffraction 
pattern and its refinement results in the Cm and R3m 
(inset) space groups at 300 K. We assumed that Fe and 
Nb are randomly distributed at the same crystallographic 
position m HU ng and allowed this sublattice to shift 
relative to Pb but did not split any sites at this stage. 
Satisfactory refinements have previously been reported 
in both R3m [7] and Cm HB; we rely on an additional 
approach to distinguish these space groups. The data are 
well explained by both structures, although the quality 
of the refinement is slightly better in Cm, but this is not 
surprising since the Cm structure model has more free 
parameters. Table [T| summarizes the refined structural 
parameters for the 300 K data in the Cm and R3m struc¬ 
tures. Results are described in the rhombohedral setting 
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of the R3m structure for comparison with Ref. [71 One 
can also convert the lattice parameters of the Cm struc¬ 
ture to pseudocubic notation to ease comparison with the 
paraelectric structure as performed in Fig. [^a) [20] . 

TABLE I. Refined structure parameters at 300 K in the Cm 
(with b as the unique monoclinic axis and cell choice 2) and 
R3m structures. The unsplit Pb position is taken as the origin 
in Cm, whereas the Fe/Nb position is set to (0.5,0.5,0.5) for 
R3m to aid comparisons to the literature. 


Cm 

ttm = 5.6819(1) A, bm = 5.6738(1) A, Cm = 


4.01202(5) A, l3 = \ 

89.896(2)° 


Ions 

X y 

z 

B (A^) 

Pb 

0.0000 0.0000 

0.0000 

1.82(8) 

Fe/Nb 0.4669(9) 0.0000 

0.5165(18) 

0.34(3) 

01 

0.4333(13) 0.0000 

0.0002(26) 

0.54(3) 

02 

0.2129(20) 0.2483(7) 0.4989(29) 

0.54(3) 

Rp- 

5.26, Rwp- 7.68, Rexp- 

4.13, x": 3.45 


R3m ttr — 4.01389(3) A, 

a = /3 = "f = 89.9223(2)° 

Pb 

0.0205(6) 0.0205(6) 0.0205(6) 

2.45(5) 

Fe/Nb 0.5000 0.5000 

0.500 

0.27(2) 

0 

0.4772(2) 0.4772(2) -0.0092(5) 

0.58(2) 

Rp. 

5.48, Rwp’- 8.01, Rexp'- 

4.13, x^: 3.76 



In the temperature-dependent data, magnetic peaks 
can be observed at d = 1.55 and 1.85 A, consistent 
with the reported G-type antiferromagnetic structure [7] . 
Several structural peaks, shown in Fig. [^b), are absent 
above the ferroelectric transition of 385 K as seen com¬ 
paring the data taken at 453 and 300 K; these peaks 
are marked by arrows. We note that we can refine the 
highest temperature data in the paraelectric phase using 
the previously reported space group Pm3m nanaEi]; 
however, as we will discuss later, there are clear signs of 
short-range structure in the total scattering data even at 
453 K. 

In order to distinguish the two candidate low- 
temperature space groups, the temperature dependence 
of the peak widths for several nuclear Bragg peaks as 
obtained from fits to a Lorentzian line shape was exam¬ 
ined, and a search was performed for peaks that are split 
in one space group but not the other. The large num¬ 
ber of reflections in this low-symmetry structure makes 
it difficult to find nonoverlapping peaks, but several good 
candidates were located. The (002) and (220) Bragg 
reflections of the Cm structure, shown in the inset of 
Fig. I^a), were found to be the most suitable: In the 
R3m space group the resulting peak contains only one 
unique reflection, (200). No clear splitting is observed at 
low temperatures, but the peak width increases markedly 
on cooling as shown in Fig. [^a), which is exactly oppo¬ 
site to what one expects from any conventional thermally 



Temperature (K) 

FIG. 2. (a) Temperature dependence of the peak width of 
the (002)/(220) Bragg peak in the Cm structure with the 
raw data at 10 K shown in the inset. Error bars are smaller 
than the symbol size. The data for the (022) peak of SrRuOa 
were taken after Ref. [22] (see the text), (b) Temperature 
dependence of Fe/Nb-0 bond lengths. The inset depicts the 
structure of PEN with the Cm space group: There are four 
different Fe/Nb-0 bonds denoted by the numbers. 


activated broadening process. Similar temperature de¬ 
pendence is observed in other peaks, which are split in 
Cm and not in R3m. As we cannot find a peak from 
a well-separated Bragg reflection in our data because of 
the low symmetry, we choose SrRuOs measured under 
almost identical conditions in order to demonstrate how 
the width of a single Bragg peak behaves as a function of 
temperature when measured at the S-HRPD beamline. 
In Fig. I^a), we plotted the (022) peak of SrRuOs [22| 
where the width of the single Bragg peak remains temper¬ 
ature independent over the measured temperature range 
as expected. The unusual peak broadening seen in the 
(002)/(220) peak of PFN is naturally explained by a split¬ 
ting of two reflections, which strengthens on cooling, and 
would appear to exclude the R3m structure. The origin 
of this unusual temperature dependence can be found 
by examining the Fe/Nb-0 bond lengths and their tem¬ 
perature dependence, shown in Fig. [^b). Two aspects 
are particularly noteworthy: First, there are already four 
different Fe/Nb-0 bond lengths even at room tempera- 












4 


ture [see the inset in Fig. [^b) for the four oxygen atoms 
with different Fe/Nb-0 bond lengths denoted by differ¬ 
ent numbers]; second, the splitting of the Fe/Nb-0 bond 
lengths only continues to grow upon cooling. Thus, the 
larger distortion of the (Fe/Nb)06 octahedron lies at the 
heart of the unusual peak broadening seen in Fig. |^a), 
favoring the Cm space group. 

The temperature dependence of the lattice constants 
and unit-cell volume are summarized in Fig. We 
also plot the monoclinic angle in pseudocubic notation. 
Square symbols represent NOVA results, whereas red cir¬ 
cles indicate those from S-HRPD. Refinements were per¬ 
formed in the Pm3m structure above Tq and the Cm 
structure below Tq. Both the lattice constants and the 
unit-cell volume follow a consistent trend without any 
clear anomaly at the antiferromagnetic transition. Re¬ 
garding the previous claim of negative thermal expan¬ 
sion m, we note that there is no evidence of such an 
effect in the data presented here, taken using either in¬ 
strument, although it is also important to note that lat¬ 
tice parameters from x-ray diffraction are normally more 
precise. It is an open question to us why our data col¬ 
lected from two different instruments do not show the 
negative thermal expansion claimed in Ref. [TTl 

In order to further analyze the temperature depen¬ 
dence of the unit-cell volume, the Debye-Griineisen for¬ 
mula [23] was used 


V{T) 


VoUjT) 
Q - hU{T) 


+ K), 


U{T) = 9NkBT 



e® - 1 


-dx. 


( 1 ) 

( 2 ) 


where U{T) is total internal energy, Q = VqBo/j, b = 
{Bq — l)/2, 7 is the Griineisen parameter, Bq and Bq 
are the bulk modulus and its first derivative with pres¬ 
sure, N represents the number of atoms in the unit cell, 
/cb is the Boltzmann constant, and Vb is volume at zero 
temperature. Fitting (shown by a solid line) results in a 
Debye temperature Qd of only 150 K, whereas b = 1.5, 
Vb = 64.582 A^, and Q = 3.87(6) x 10“^^ J are simi¬ 
lar to parameters used for other materials [24ll26] . We 
note that the resulting Debye temperature is abnormally 
small, leading to a bulk modulus around 600 GPa for a 
typical Griineisen parameter of order one, which is high 
for an oxide: For example, we earlier found bulk mod¬ 
uli of 120 GPa for YMnOs [27| while 250 GPa has been 
reported for MgSiOs [28]. Note that this assumption 
about the Griineisen parameter is crude, and so our es¬ 
timate of the bulk modulus should be taken with cau¬ 
tion. Moreover a true test of the bulk modulus needs 
to be performed by measuring the volume measurement 
under pressure. On the other hand, fixing the Debye tem¬ 
perature to the more typical 430 K led to a noticeably 
poorer fit (dashed line) with b = 1.5, Vb = 64.595 A^, and 
Q = 2.9(1) X 10“^^ J. This could be consistent with an 


anomaly below Tn but opposite in sign to that reported 
previously. 

Figure |^d) shows the temperature dependence of 
atomic displacement parameters for each atom. The 
thermal parameter of Pb is much larger than those of 
O and Fe/Nb, almost comparable to those of other Pb- 
containing ferroelectric materials m and consistent with 
a previous report El- The temperature dependence was 
modeled using 



where i represents each atomic species (Pb, Fe/Nb, O), 
and A* = 36.39/M*0^ is related to zero-point energy 
of the atoms concerned with atomic mass M* m- How¬ 
ever, it was necessary to add constant offsets to A to 
achieve better agreement. The final result, denoted by 
lines in Fig.j^d), indicated effective Debye temperatures 
of 150 K for Pb, 680 K for O, and 370 K for Fe/Nb. 

Using the structural information, we then searched for 
possible experimental evidence of magnetoelectric cou¬ 
pling. PFN is known to exhibit a frequency-dependent 
anomaly in the dielectric constant at Tn, but little else is 
known of the magnetoelectric coupling. From the high- 
accuracy structural information presented here, it is pos¬ 
sible to calculate the expected ferroelectric polarization, 
starting from the centrosymmetric Pm3m space group. 
Using a similar approach, we were previously able to show 
a negative magnetoelectric coupling in one of the best- 
studied multiferroic compounds, BiFeOs [29]. For sim¬ 
plicity, the nominal valences of the atoms were used as a 
starting point: Pb^+, Fe^+, Nb^+, and 0^“. Since the 
actual valences, so-called Born effective charges, which 
include dynamic terms, might well differ from the nom¬ 
inal charge values [30], the discussion below may not be 
fully quantitative, but it should be qualitatively correct. 
Using the above assumptions, the electric polarization 
was calculated with respect to the paraelectric Pm3m 
phase. As shown in Fig. |^e), the estimated polariza¬ 
tion at room temperature is around 20/iC/cm^, which 
is twice that reported EE] in poly crystalline samples. 
Considering the rough assumptions made and the fact 
that the experiments were performed on powder sam¬ 
ples, we consider this degree of deviation acceptable. 
The temperature dependence was then parametrized us¬ 
ing Ginzburg-Landau analysis, assuming that the electric 
polarization follows the usual first-order temperature de¬ 
pendence: F{P) = + ^P^ + ^P^, where F is 

the Ginzburg-Landau free energy and P is the polariza¬ 
tion. A fit to all points produced an unphysically high 
Curie temperature, whereas constraining the Curie tem¬ 
perature led to clear systematic trends in the residuals; 
a fit to only the points above the Neel transition, shown 
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FIG. 3. Temperature dependence of (a) lattice parameters, (b) monoclinic angle, (c) unit-cell volume, (d) thermal parameter, 
and (e) calculated polarization. The refined unit-cell volume (symbols) is shown together with the theoretical curve obtained 
from the Debye-Griineisen formula, and polarization was calculated from the rehned atomic positions. 


in Fig. [3];e), predicted a Curie temperature of 380 K in 
close agreement with experiment, but the data appear to 
deviate from this fit below Tn- This discrepancy, marked 
by shading in the figure, is the clearest experimental ev¬ 
idence yet of the magnetoelectric coupling in PFN. 

For further analysis of the short-range structure, 
Fig. Qa) shows how the radial distribution function G(r) 
calculated from the NOVA data changes gradually with 
temperature. Of note are the several regions where 
peaks’ intensities increase on cooling. At these short 
length scales, corresponding to large momentum transfers 
Q, we can ignore contributions from magnetic scattering 
because they are usually small compared to nuclear scat¬ 
tering m- As shown in Fig. [^b), the long-range struc¬ 
ture up to 100 A can be well explained by the Cm space 
group. 

First, we consider the short-range structure at high 
temperatures above the ferroelectric transition. Using 
the data taken at 453 K, we tested four different mod¬ 
els to explain the data and clarify the local structure, all 
variants of the Cm structure, and all explained in more 
detail below: (A) one without local short-range struc¬ 
ture, (B) one with Pb displacements as considered pre¬ 
viously [14], (C) one with Pb disorder, and (D) one with 
both Pb displacement and O disorder. For all models, Fe 
and Nb are treated as an average pseudoatom as in the 
refinement. Models B-D take model A as their starting 


point and allow additional freedom in the atomic posi¬ 
tions. In departing from model A, C centering was not 
retained. As shown in Fig. Qc), the agreement is indeed 
improved by including Pb displacement and O disorder 
in our model calculation. Upon closer inspection, model 
A is found not to satisfactorily explain the intensities of 
peaks at r = 2 and 3.5 A, and agreement with experimen¬ 
tal results is poor at longer distances. It is particularly 
surprising that, even in the nominally paraelectric phase 
well above the ferroelectric transition, there is clear ev¬ 
idence of short-range structure. We note that similar 
behavior is often found in relaxor ferroelectric materi¬ 
als [32l[33]. Thus, our result adds further weight to the 
previous claim that PFN is a rare relaxor multiferroic. 

Having concluded that the simplistic model A does 
not adequately explain even the high-temperature local 
structure, we now consider other models which might de¬ 
scribe the data at all temperatures. Given lead’s large 
thermal parameter even after accounting for the offset 
between the Pb and Fe/Nb sublattices in the Cm struc¬ 
ture (model A) as shown in Fig. [^d), it is natural to 
think of further Pb displacements as an origin of short- 
range structure as is often the case in closely related 
Pb-containing perovskite systems [T21|5i] and as consid¬ 
ered previously [T4|. At 62 K, model B requires that the 
lead atoms shift by (—0.019,0.03,0.012) in fractions of the 
pseudocubic unit cell, less than previously reported mi. 











6 



FIG. 4. (a) Radial distribution function G{r) at several temperatures, calculated from NOVA results. Gradual changes may 
be observed with temperature at several regions of interest (see the text), (b) Radial distribution function G{r) is compared 
with fitting results over a wide length-scale range for the 62 K data. The symbols represent the data points, whereas the line 
is the fitting results with the difference curve shown at the bottom. The local structure is compared with several models at (c) 
453 K and (d) 62 K. Gircles are calculated from experimental data, and lines are fits to several models: Model A assumes no 
local structure; model B assumes Pb displacements; model G considers Pb disorder; and model D combines Pb displacements 
with O disorder. Traces have been offset vertically for clarity. 


Moving the Pb atoms introduces several peaks between 
r = 3 and 4 A, but their intensities are much smaller than 
in the experimental results, and the peak at r = 2.5 A 
is not reproduced. Above the ferroelectric transition, 
symmetry constrains the Pb displacement to be ( 77 , 77,77), 
where 77 = 0.02 at 453 K. Displacing the Pb atoms in con¬ 
cert, which preserves the unit cell, is quite a simplistic 
model, and allowing the Pb atoms more freedom of move¬ 
ment should better model the real system. Accordingly, 
model C was introduced with the possibility of random 
shifts in the Pb positions. 

Model C simulates the effect of Pb disorder using a 
4x4x4 supercell in which each lead atom was allowed 
to move freely from its original position. Although this 
model was able to better model the experimental radial 
distribution function, at least at low temperatures, the 
average Pb displacement from its already displaced Cm 
position is 0.2 A in this model. This is very large com¬ 
pared to Pb’s atomic thermal motion in Fig. |^d), ap¬ 
proaches the upper limit observed in the strongest Pb- 
based ferroelectrics and BiFeOs with one of the high¬ 
est polarization values of r\j 86 /iC/cm^ [29], and would 
suggest an electric polarization far exceeding that re¬ 
ported mill- The comparative success of model C rela¬ 
tive to model B is most likely a result of the significantly 
higher number of free parameters. Although model C is 
ruled out on the basis of the extremely large displace¬ 
ments compared to the reported polarization values, the 
fact that it required unphysical displacements of Pb im¬ 


plies that other atoms must also be considered. It has 
already been shown that the Fe/Nb atoms do not ex¬ 
hibit any significant displacement from their ideal posi¬ 
tions [H]. That leaves oxygen, which has a larger ther¬ 
mal parameter and which is well known to depart from 
its ideal position in most perovskite structures, includ¬ 
ing in several Pb-containing phases which exhibit oxygen 
displacement [35H37] or the random rotation of oxygen 
polyhedra [33]. Model D allows the oxygen atoms to 
shift. 

A further justification for considering O disorder comes 
from the reported random distribution of neighboring 
Fe3+ and Nb^+ and the lack of any features in our 
data that would suggest even short-range cation order. 
The magnetic moments calculated from the long-range 
refinement are 3.4(6) fiB at 62 K (NOVA) and 3.3(3) jiB 
at 10 K (S-HRPD), corresponding to the high spin state 
of Fe^+, which would make the two cations’ sizes essen¬ 
tially identical and eliminate the common source of dis¬ 
order. The next most probable reason for oxygen atoms 
to depart from their ideal positions is the displacements 
of the Pb^+ cations from their already-shifted positions, 
which would be expected to cause small rotations of the 
(Fe,Nb)06 octahedra. This would be frustrated, unless 
the Pb displacements ordered in one of a few very spe¬ 
cific ways which would be expected to produce additional 
nuclear reflections. An orthogonal possibility arises from 
the electric-field gradients created by having a random 
distribution of Fe^+ and Nb^+ cations. An oxygen atom 
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situated between an iron and a niobium site feels an 
electric-field gradient that will polarize it. Meanwhile, 
Nb^+, having the electronic configuration of krypton, will 
be unable to share electron density with oxygen, which 
3d^ Fe3+ can, using its weakly antibonding e* orbitals. 
In fact, iron, being a 3d transition metal, should have 
significant on-site electron repulsion and may benefit by 
donating some of its electron density to its oxygen lig¬ 
ands. All of these electronic effects would produce a net 
flow of electron density away from iron sites and toward 
niobium sites with many oxygen atoms caught in the mid¬ 
dle. In model D, we combined lead displacement (model 
B) with oxygen disorder using a 4 x 4 x 4 supercell. Here 
we allowed oxygen atoms to move freely while all lead 
atoms were constrained to move together. As one can 
see in Figs. [^c) and Bd), this model apparently fits the 
experimental results very well. The fitted atomic dis¬ 
placement in fractions of the pseudocubic unit cell for 
lead is (0.010,-0.004,0.008) at 62 K and (—0.0051) at 
453 K, which is reduced compared to model B. For oxy¬ 
gen, the average shifts relative to the original site are 
(0.0195) at 62 K and (0.0203) at 453 K in fractional co¬ 
ordinates, which is around 0.08 A. As is clear in Fig. 
of the models considered here, model D is the most suc¬ 
cessful at describing the experimental results, and it does 
so in a physically reasonable manner. 

In Fig. we deconvoluted the contributions of corre¬ 
lations from each pair of atoms using models A and D for 
the 62 and 453 K data shown in Fig. In each panel, 
the top six lines show correlation functions between in¬ 
dividual pairs of atoms as labeled. Below this are the 
experimental results and best fit and finally the residu¬ 
als. A comparison of panels (a) with (b) in Fig. [^reveals 
that 0-0 correlations develop a great deal of additional 
structure at both temperatures, explaining well the ad¬ 
ditional peaks at r = 2.5-4 A. 

CONCLUSION 

In conclusion, we have undertaken both high-resolution 
powder-diffraction and total scattering experiments on 
PFN in order to shed light on the unusual relaxor mul- 
tiferroicity in this material, focusing on four key aspects 
of the physical properties: (i) determination of the low- 
temperature structure, (ii) the experimental examination 
of the reported negative thermal expansion behavior, (iii) 
possible structural evidence of magnetoelectric coupling, 
and (iv) the nature of the relaxor behavior. By com¬ 
bining both long- and short-range diffraction data, it is 
possible to provide new insight on each of the four is¬ 
sues. First, the high-resolution neutron-diffraction data 
are more consistent with the Cm space group at low tem¬ 
peratures as reported based on synchrotron data m and 
exclude the previously proposed R3m [7] structure. How¬ 
ever, there are no signs of the negative thermal expansion 


(a) 453 K-A (b) 453 K-D 




FIG. 5. Atomic position correlations at 453(62) K with (a) 
[(c)] no local structure and (b) [(d)] both Pb displacement and 
O disorder. The top six lines show atom-atom radial correla¬ 
tion functions as labeled, followed by the experimental results 
(circles), best fit (line), and finally its residuals (bottom line). 
Traces have been offset vertically for clarity. 

reported in Ref. [TT] with the data instead being well fit 
by the conventional Debye-Griineisen formula. We suc¬ 
ceeded in finding a structural signature of the magne¬ 
toelectric coupling by Ginzburg-Landau analysis of the 
calculated electric polarization. Finally, the data clearly 
show that over the entire temperature range covered in 
this study there exists a short-range structure in addi¬ 
tion to the global structure of Pm3m or Cm. This lo¬ 
cal structure occurs through Pb displacement as well as 
O disorder. That this local structure is present in PFN, 
even in the paraelectric phase, as seen in relaxor ferroelec¬ 
tric materials may hold the key to the strong frequency 
dependence seen in the dielectric constant of PFN. 
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